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• Premise of the study: We isolated 15 polymorphic microsatellite markers from Vachellia farnesiana for use in population 
genetic studies to determine the native range of the species. 

• Methods and Results: Initially, 454 shotgun sequencing was used to identify and design primers for 68 microsatellite loci. 
Of these, we trialed 47 loci in the target species, and 42 (89%) amplified a product of expected size. Fifteen of the 47 loci were 
screened for variation in 21 individuals from the native range of V. farnesiana in southern Mexico and 20 from northwestern 
Australia. Fourteen loci were polymorphic, with observed heterozygosity ranging from 0.026 to 1.00 (mean = 0.515) and two 
to 12 alleles per locus (average = 5.2). Cross-amplification was successful in four to 11 loci in three other Vachellia species. 

• Conclusions: The new microsatellite loci will be useful in understanding genetic variation and investigating the role of human- 
mediated dispersal in the current distribution of V. farnesiana. 

Key words: 454 GS-FLX; Acacia farnesiana; cross-species transferability; microsatellites; shotgun sequencing; Vachellia 
farnesiana. 



Mimosa bush, Vachellia farnesiana (L.) Wight & Arn. (syn- 
onym Acacia farnesiana (L.) WillcL), is a woody mimosoid le- 
gume with a pantropical distribution. It has several common 
names in its native range, including mimosa bush, sweet acacia, 
cassie, and huizache. Acacia Mill., if treated in the broad sense 
(sensu lato [s.l.]), is a large polyphyletic genus, with at least five 
lineages that may be recognized as genera: Acacia sensu stricto 
(s.s.), Acaciella Britton & Rose, Mariosousa Seigler & Ebinger, 
Senegalia Raf., and Vachellia Wight & Arn. (Maslin, 2008; 
Bouchenak-Khelladi et al., 2010). The genus Vachellia is com- 
posed of a predominantly African clade and a predominantly 
American clade (Bouchenak-Khelladi et al., 2010). Vachellia 
farnesiana is part of the American clade, but has a distribution 
that extends well beyond the Americas, and it is considered in- 
vasive in some countries. Its arrival date in Australia, and hence 
its status as native or alien, remains unknown, but V. farnesiana 
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may have arrived prior to European colonization (Bean, 2007). 
The Spanish and Portuguese introduced the species to Europe 
in the 17th century. At this time, the two countries had a strong 
colonial presence around the Indian Ocean, through which fur- 
ther dispersal of the plant was possible. However, natural ocean 
currents and pre-European indigenous traders may have played 
a role in earlier dispersals. Genetic data from V. farnesiana may 
be useful in determining the dispersal pathways of this plant to 
populations outside of the Americas. Microsatellite markers 
have been developed previously for the invasive V. nilotica (L.) 
P. J. H. Hurter & Mabb. (Wardill et al, 2004). However, only a 
total of five loci were developed, and it is unknown how many 
of these will cross amplify in V. farnesiana. It was, therefore, 
necessary to develop new markers for V. farnesiana to facilitate 
our investigations of population genetics and plant dispersal out 
of the native range. 

METHODS AND RESULTS 

Genomic DNA (5 jig) was isolated from one individual of V. farnesiana 
from silica gel-dried leaves with the QIAGEN DNeasy Plant Mini Kit (QIA- 
GEN, Valencia, California, USA) as per the manufacturer's protocol. The DNA 
was sent to the Australian Genomic Research Facility (AGRF) in Brisbane, 
Australia, for shotgun sequencing on a Titanium GS-FLX (454 Life Sciences, a 
Roche Company, Branford, Connecticut, USA) following Gardner et al. (201 1). 
The sample occupied 12.5% of a plate and produced 59,289 individual se- 
quences, with an average fragment size of 307 bp; 1.9% of the sequences con- 
tained microsatellites. The raw data from shotgun sequencing were deposited in 
the Dryad Digital Repository (doi:10.5061/dryad.jdl83; Meglecz et al., 2012). 
We used the program QDD version 1.3 (Meglecz et al., 2010) to screen the raw 
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sequences with eight or more di-, tri-, tetra-, or pentabase repeats. Redundant 
sequences were removed and primers were designed with a specified PCR 
product length of 80-480 bp using Primer3 (Rozen and Skaletsky, 2000) in 
QDD; default settings were maintained for all parameters except product length. 
The software identified and designed primers for a total of 68 loci, of which 47 
contained simple repeats and 21 contained tandem repeats (Table 1). 

We followed the procedure outlined in Gardner et al. (201 1) for further de- 
velopment of the 47 loci containing simple repeats. The 47 loci were trialed for 
amplification using seven V. farnesiana individuals, each from a different pop- 
ulation (Appendix 1), and 10-jiL reactions containing lx buffer, 0.5 U HotStar- 
Taq DNA polymerase (QIAGEN), 1.5 mM MgCl 2 , 0.25 mM of each dNTP, 
250 nM each forward and reverse locus-specific primer, and 10-50 ng genomic 
DNA. The following PCR conditions were used: 95°C for 15 min; followed by 
28 cycles at 95°C for 30 s, 58°C for 90 s, and 72°C for 30 s; and a final elonga- 
tion step at 60°C for 30 min. PCR products were visualized on a 1.5% agarose 
gel stained with ethidium bromide. Twenty-eight loci amplified a product of the 
expected size for all seven samples, with no unexpected secondary bands. 
These 28 loci were tested for polymorphism using forward primers tagged with 
454A sequence tags and 454A sequencing tags labeled with either 6-FAM, 
NED, HEX, or PET (Applied Biosystems, Foster City, California, USA) fol- 
lowing the method of James et al. (2011) and were run by Macrogen (Seoul, 
Korea) on a 3730x1 DNA sequencer (Applied Biosystems) with a GeneScan 
500 LIZ Size Standard (Applied Biosystems). Of the 28 loci tested, 26 loci 
(93%) were polymorphic, one (3.6%) was monomorphic, and one (3.6%) did 
not amplify for all samples under these conditions. Of the 26 polymorphic loci, 
1 1 (42%) produced alleles that were affected by stuttering or amplified weakly 
and were removed from further consideration. The remaining 15 (54%) poly- 
morphic loci (Table 1) were screened for variation in 20 recently collected in- 
dividuals from a single population from southern Mexico, one herbarium 
specimen also from southern Mexico, and 20 recently collected individuals 
from northwestern Australia (Table 2), with DNA isolation, PCR, and fragment 
length analysis as described above. For each locus, we calculated the number 
and range of alleles, observed (H 0 ) and expected heterozygosity (H e ), and de- 
viation from Hardy-Weinberg equilibrium (HWE) using GenAlEx (Peakall 



and Smouse, 2006). P values from HWE tests were adjusted for multiple tests 
of significance using the sequential Bonferroni method (Holm, 1979). The 
number of alleles per locus ranged from one to 12 across these 41 individuals, 
and H e ranged from 0 to 0.84. Within the Mexican samples, seven polymorphic 
loci were in HWE, five significantly deviated from HWE, and three were mono- 
morphic. Within the Australian samples, nine polymorphic loci were in HWE, 
five significantly deviated from HWE, and one was monomorphic (Table 2). 
We used MICRO-CHECKER 2.2.3 (van Oosterhout et al., 2004) to check each 
locus for further evidence of null alleles, scoring error due to stuttering, and 
large allele dropout. Four loci (Af03, Af47, Af32, Af26) showed significant 
null allele frequencies at the target site, or evidence of scoring error due to stut- 
tering. None of the loci showed evidence of large allele dropout. We checked 
all pairs of loci for linkage disequilibrium in GENEPOP and none were signifi- 
cant after sequential Bonferroni adjustment. 

Primers for the 15 selected loci were also tested for amplification and cross- 
species transferability in 12 individuals of V. nilotica (8 recently collected and 
4 herbarium specimens), two herbarium specimens of V. aroma (Gillies ex 
Hook. & Arn.) Seigler & Ebinger, and one each of the Australian species 
V. ditricha (Pedley) Kodela and V. suberosa (A. Cunn. ex Benth.) Kodela 
(Appendix 2). Isolation of DNA, PCR, and fragment analysis were as described 
above. Thirteen of the 15 loci amplified successfully in the majority of indi- 
viduals of V. nilotica, and eight of these were polymorphic for the small number 
of individuals examined. Amplification success was lower for the remaining 
species (5-12 of 15 loci), possibly due to the use of DNA isolated from her- 
barium specimens. 



CONCLUSIONS 

These markers will be used to document the genetic diversity 
of V. farnesiana and to investigate the dispersal pathways lead- 
ing to its current pantropical distribution. Given the successful 



Table 1. Characterization of 15 polymorphic micro satellite loci of Vachellia farnesiana. 



Locus a Primer sequences (5'-3') GenBank accession no. Repeat motif Allele size (bp) T & (°C) 



Afl8 v 


F 
R 


GCCACAACTAAAGTCATATCACCA 
CCTTCTTACGCTCCATGATTC 


KF030919 


(TA) 9 


108 


58 


Af24 p 


F 
R 


CATGGCCTATTTCCACCACT 
TTGGTGCAATTGATAGCGTT 


KF030921 


(AT) 9 


94 


58 


Af05 p 


F 
R 


TTGGACATTCCAATTGAGATTATTA 
AGCAGGAACTTGCTTAGATGC 


KF030916 


(TG) 8 


118 


58 


Af38 F 


F 
R 


GATTGCTATGTCATCTCCCTCC 
GTGCGAGATCTATCGACGAC 


KF030926 


(GT) 10 


98 


58 


Afl9 F 


F 
R 


ACTTCGAGATGAACCTCCCA 
CGAGACCCAAATCAGTCGAT 


KF030920 


(AT) n 


106 


58 


Af32 N 


F 
R 


CAGTTCAAACTATCATCTCTATTCACA 
GTGATATGTTTACGGTGCCGA 


KF030925 


(AT) 8 


90 


58 


Af25 N 


F 
R 


GATGGCGGCAACACAGTAT 
AAGT GAAC AAT AT T GAAGCGCA 


KF030922 


(CT)io 


109 


58 


Af03 N 


F 
R 


TTAATGCAATTGGGAATCACTT 
GACACTCCCACCTGTATCGG 


KF030915 


(GA) 15 


150 


58 


Afl0 F 


F 
R 


GAAGTTATTCTTAATTGCTACCATTCC 
TTGACCAACTCTACTCTTAATTGATTG 


KF030917 


(AC) 12 


91 


58 


Af26 F 


F 
R 


CAGCTCGATAGCTAAACAAGGA 
GGTGTTTGGATGGAAGTTCG 


KF030923 


(CA) 10 


108 


58 


Af47 F 


F 

R 


CCTGAGACAGTTGTGTTTGATTG 
ATCATGCCTTGTCAGCATCC 


KF030929 


(AQn 


121 


58 


Afl4 N 


F 

R 


ATTACACCACTCGGTCGGTC 
CCCATCTTCTCCAGCATCAT 


KF030918 


(AAG) 5 


90 


58 


Af29 N 


F 

R 


GGAATCCAATGTATTTGGCG 
AGGTTCACAAGGCAACCTGT 


KF030924 


(AT) 8 


109 


58 


Af42 N 


F 

R 


AAAC T C AAT AAC T T GC T T AAC T GAAA 
CCAATTTGCTTGCTTGACTTG 


KF030927 


(TC) 5 


120 


58 


Af46 N 


F 

R 


TGAAGAATAATAGCTAGCGGCTG 
TGAGAAGGCCCAATGAAATC 


KF030928 


(AG) 9 


91 


58 



Note: T a = annealing temperature. 

a Superscripts F, N, V, and P indicate loci were 5' labeled with the dyes 6-FAM, NED, VIC, and PET, respectively. 
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cross-amplification of these loci for a broad range of Vachellia 
species, the primers may be useful for studies of the genetic 
diversity of other Vachellia species. 
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Appendix 1. Locality data for the seven individuals of Vachellia farnesiana used in the initial screening of 47 loci. 



Country 


Collection locality 


Geographic coordinates 


USA 


Arizona: Maricopa County 


33.2932°N, 112.428°W 


Mexico 


Veracruz: Los Negritas 


18.8383°N, 96.07°W 


Mexico 


San Luis Potosf 


22.2°N, 101°W 


Madagascar 


Antsiranana: Diana: south of Diego Suarez 


12.4321°S, 49.3567°E 


Madagascar 


Nosy Be 


13.3833°S,48.2°E 


Australia 


Queensland: 3 1 km W of Cloncurry 


20.7584°S, 140.2327°E 


Australia 


Western Australia: 180 km E of Halls Creek 


17.944°S, 128.8816°E 



Appendix 2. Voucher information for Vachellia species used in this study. 



Species 


Voucher specimen accession no. a 


Collection locality 


Geographic coordinates No. of individuals 


Vachellia farnesiana 


MO6178804 b 


Oaxaca, Mexico 


16.302°N, 96.286°W 1 


V farnesiana 


K.L. Bell 128 c 


Katherine, Northern Territory, Australia 


14.463°S, 132.259°E 1 


V farnesiana 


ASU 279693 d 


Maricopa County, Arizona, USA 


33.2932°N, 112.428°W 1 


V farnesiana 


MEL 2370354A 


Diana, Antsiranana, Madagascar 


12.4321°S, 49.3567°E 1 


V aroma 


MEL 2263911 


Bolivia 


20.105°S, 63.487°W 1 


V aroma 


MEL 2263912 


Bolivia 


17.9°S, 64.558°W 1 


V ditricha 


MEL 2066644 


Wyndham-East Kimberley, Western Australia, Australia 


16.3839°S, 126.4975°E 1 


V. suberosa 


MEL 2066645 


Wyndham-East Kimberley, Western Australia, Australia 


16.3839°S, 126.4975°E 1 


V nilotica 


MEL 260774 


Queensland, Australia 


21.267°S, 141.3°E 1 


V nilotica 


MEL 2080462 


Queensland, Australia 


20.05°S, 148.25°E 1 


V nilotica 


MEL 2204859 


Western Australia, Australia 


15.803°S, 128.75°E 1 


V nilotica 


MEL 2293312 


Queensland, Australia 


23.446°S, 150.439°E 1 


a Lodged at the National Herbarium of Victoria (MEL), except where noted. 
b Lodged at the Missouri Botanical Garden (MO). 
c Lodged at MEL, but not yet accessioned. 



d Lodged at the Arizona State University Vascular Plant Herbarium (ASU). 
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